A First Measurement of Low x Low Q 2 Structure Functions in Neutrino Scattering 
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Neutrino structure function measurements in the low 
Bjorken x, low Q 2 region can be used to study the axial- 
vector component of the weak interaction as well as to 
test the limits of parton distribution universality. In this 
paper, we present a first measurement of the structure 
function F2 in neutrino scattering, from the CCFR data, 
for Q 2 < 1 GeV 2 and 0.0045 < x < 0.035. In this region 
where perturbative and non-perturbative QCD meet, we 
present a parameterization of the data which allows us to 
test the partially conserved axial current (PCAC) limit 
of F2 in neutrino scattering. 

The universality of parton distributions can be tested 
by comparing neutrino scattering data to charged lepton 
scattering data. Past measurements for 0.0075 < x < 0.1 
and Q 2 > 1.0 GeV 2 have indicated that F% differs from 
F£ by 10-15% ||. This discrepancy has been partially 
resolved by recent analyses of F% at Q 2 > 1.0 GeV 2 
While we expect and have now observed that parton dis- 
tribution universality holds in this region, this need not 
be the case at lower values of Q 2 . Deviations from this 
universality at lower Q 2 are expected due to differences in 
vector and axial components of electromagnetic and weak 
interactions. In particular, the electromagnetic interac- 
tion has only a vector component while the weak interac- 
tion has both vector and axial- vector components. Vector 
currents are conserved (CVC) but axial-vector currents 
are only partially conserved (PCAC). Adler f| proposed 
a test of the PCAC hypothesis using high energy neu- 
trino interactions, a consequence of which is the predic- 
tion that F2 approaches a non-zero constant as Q 2 — > 
due to U(l) gauge invariance. A determination of this 
constant is performed here by fitting the low Q 2 data to 
a phcnomenological curve developed by Donnachie and 
Landshoff 0. 

The differential cross-sections for the charged- 
current process (v^) + N — > (f~ (/^ + ) + X in the limit 
of negligible quark masses and neglecting lepton masses, 
in terms of the Lorentz- invariant structure functions F2, 
2xF\, and xF$ are: 



da u ' v _ G 2 F ME y 
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where Gf is the weak Fermi coupling constant, M is 
the nucleon mass, E v is the incident v energy, Q 2 is the 
square of the four-momentum transfer to the nucleon, 
the scaling variable y — Eh ad /E v is the fractional en- 
ergy transferred to the hadronic vertex with Eh ad equal 
to the measured hadronic energy, and x — Q 2 /2ME v y, 
the Bjorken scaling variable, is the fractional momen- 
tum carried by the struck quark. The structure func- 
tion 2xF\ is expressed in terms of F2 by 2xFi(x,Q 2 ) = 

F 2 (x,Q 2 ) x ^X+RixG?) > where R = ft is the ratio 
of the cross-sections of longitudinally to transversely- 



polarized W bosons. In the leading order (LO) quark- 
parton model, F2 is the sum of the momentum densities 
of all interacting quark constituents, and xF 3 is the dif- 
ference of these, the valence quark momentum density; 
these relations are modified by higher-order QCD correc- 
tions. 

The v deep inelastic scattering (DIS) data were col- 
lected in two high-energy high-statistics runs, FNAL 
E744 and E770, in the Fermilab Tevatron fixed-target 
quadrupole triplet beam (QTB) line by the CCFR collab- 
oration. The detector |8||| consists of a target calorime- 
ter instrumented with both scintillators and drift cham- 
bers for measuring the energy of the hadron shower, 
Eh ad, and the muon angle, 9^, followed by a toroid spec- 
trometer for measuring the muon momentum p^. There 
are 1,030,000 events and 179,000 F M events in the 
data sample after fiducial volume and geometric cuts, 
and kinematic cuts of p M > 15 GeV, 6^ < 150 mr, Eh ad 
> 10 GeV, and 30 < E v < 360 GeV. These cuts were 
applied to select regions of high efficiency and small sys- 
tematic errors in reconstruction. 

The structure function F2 in Eq. (|l|) can be calculated 
from the observed number of and events combined 
with the V/j, and fluxes. The ratio of fluxes between dif- 
ferent energies in v mode and that between V and v mode 
was determined using the events with Eh ad < 20GeV 
]To| |l2"[ . The overall normalization of the flux was con- 
strained such that the measured total neutrino-nucleon 
cross-section for vs equaled the world average cross- 
section for isoscalar-corrected iron target experiments, 
a vFe jE = (0.677 ±0.014) x 10~ 38 cm 2 /GeV §H$L§ and 
for Fs equaled the world average cross-section including 
this experiment for isoscalar-corrected iron target experi- 
ments, <j VFe /E = (0.340 ±0.007) x 10" 38 cm 2 /GeV. Neg- 
ligible corrections for non-isoscalarity of the iron target 
and the mass of the W boson propagator are applied. 

Sources of systematic error on F2 arise from limitations 
of the models used for corrections and from the level of 
our knowledge of the detector calibration. Muon and 
hadron energy calibrations were determined from test 
beam data collected during the course of the experiment 
p,Pj. For acceptance, smearing, and radiative correc- 
tions we chose an appropriate model for the low x, low 
Q 2 region, the GRV jyj model of the parton distribution 
functions. The GRV model is used up to Q 2 = 1.35 GeV 2 
where it is normalized to a LO parameterization [[15[ used 
above this. Inclusion of the GRV model in the radiative 
correction calculation caused a systematic decrease in F2 
by as much as 10% in the lowest x bin, decreasing to 
1-2% at x = 0.015 as compared to the effects of the LO 
model used in the previous analysis jl^] . Due to the sys- 
tematic uncertainty in the model at low x, the radiative 
correction error is 3% in the lowest x bin. A correction is 
applied for the difference between xF£ and xFg, deter- 
mined using a LO calculation of AxF$ = xFg— xF%. The 
recent CCFR AxF% measurement [El is higher than this 
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LO model [15| and all other current LO and NLO theoret- 
ical predictions in this kinematic region. An appropriate 
systematic error is applied to account for the differences 
between the theory and this measurement. Finally, a 
systematic error is applied to account for the uncertainty 
in the value of R, which comes from a global fit to the 
world's measurements [ fl6| . 

In previous analyses a slow rescaling correction was 
applied to account for massive charm effects. This is 
not applied here since the corrections are model depen- 
dent and uncertain in this kinematic range. As a re- 
sult, neutrino and charged lepton DIS data must be com- 
pared within the framework of charm production mod- 
els, accomplished by plotting the ratio of data to theo- 
retical model. The theoretical calculation corresponding 
to the CCFR data employs NLO QCD including heavy 
flavor effects as implemented in the TR-VFS(MRST99) 
scheme [[l7]jll|. The theoretical calculation correspond- 
ing to NMC and E665 data is determined using TR- 
VFS(MRST99) for charged lepton scattering. Other 
theoretical predictions such as ACOT-VFS(CTEQ4HQ) 
HQ and FFS(GRV94) [§l| do not significantly change 
the comparison. 

The combination of the inclusion of the GRV model at 
low x and low Q 2 , its effect on the radiative corrections, 
and removal of the slow rescaling correction help to re- 
solve the longstanding discrepancy between the neutrino 
and charged lepton DIS data above x = 0.015. F2 is 
plotted in Figure [l]. Errors are statistical and systematic 
added in quadrature. A line is drawn at Q 2 — 1 GeV 2 to 
highlight the kinematic region this analysis accesses. Fig- 
ure U compares Fi (data/theoretical model) for CCFR, 
NMC, and E665. There is agreement to within 5% down 
to x = 0.0125. Below this, as x decreases, CCFR F2 
(data/theory) becomes systematically higher than NMC 
F2 (data/theory). Differences between scattering via the 
weak interaction and via the electromagnetic interaction 
as Q 2 — > may account for the disagreement in this re- 
gion. 

In charged lepton DIS, the structure function F2 is 
constrained by gauge invariance to vanish with Q 2 as 
Q 2 — » 0. Donnachie and Landshoff predict that in the 
low Q region, Fg will follow the form §: 

Q 2 
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However, in the case of neutrino DIS, the axial compo- 
nent of the weak interaction may contribute a nonzero 
component to F2 as Q 2 approaches zero. Donnachie and 
Landshoff predict that F% should follow a form with a 
non-zero contribution at Q 2 = 0: 
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rameterization of SLAC Fe/D data JTq|, we do a com- 
bined fit to the form predicted for /i DIS and extract the 
parameter A = 0.81 ± 0.02 with X 2 /DOF = 27/17. Re- 
sults of fits in each x bin for each experiment are shown 
in Table | for comparison to parameters in the CCFR 
fit. The error on A is incorporated in the systematic er- 
ror on the final fit. Inserting this value for A into the 
form predicted for i/N DIS, we fit CCFR data to extract 
parameters B, C, and D, and determine the value of F2 
at Q 2 = 0. Only data below Q 2 = 1.4 GeV 2 are used 
in the fits. The CCFR ir-bins that contain enough data 
to produce a good fit in this Q 2 region are x = 0.0045, 
x = 0.0080, x = 0.0125, and x = 0.0175. Figure | and 
Table || show the results of the fits. Error bars con- 
sist of statistical and systematic terms added in quadra- 
ture but exclude an overall correlated normalization un- 
certainty of 1-2%. The values of F 2 at Q 2 = GeV 2 
in the three highest x-bins are statistically significant 
and are within 1 a of each other. The lowest x bin has 
large error bars but is within 1.5cr of the others. Taking 
a weighted average of the parameters B,C,D, and F2 
yields B = 1.53 ± 0.02, C = 2.31 ± 0.03, D = 0.48 ± 0.03, 
andF 2 (Q 2 = 0) = 0.21±0.02. Figure | shows F 2 (Q 2 = 0) 
for the different x bins. Inclusion of an x dependence of 
the form x 13 does not change the overall fits or x 2 s. How- 
ever, unlike in charged lepton scattering, the Donnachie 
and Landshoff mass parameter, B, appears to depend on 
x, with higher values corresponding to higher x. Thus, 
F2 at higher x approaches F2(Q 2 = 0) more slowly than 
at lower x. 

In summary, a comparison of F2 from neutrino DIS 
to that from charged lepton DIS shows good agreement 
above x = 0.0125, but shows differences at smaller x. 
This low x discrepancy can be explained by the different 
behavior of F2 from v DIS to that from e//j, DIS as Q 2 — ► 
0. CCFR F% data favors a non-zero value for F2 as Q 2 — > 
0. 
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FIG. 2. F 2 data/theory from CCFR v-Ye DIS compared to 
F2 from NMC and E665 DIS. Errors bars are statistical and 
systematic added in quadrature. Theoretical predictions are 
those of TR-VFS(MRST99). 
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FIG. 1. CCFR F 2 at low x, low Q 2 . Data to the left of the 
vertical line at Q 2 = 1.0 represent the new kinematic regime 
for this analysis. 



TABLE I. Results for NMC and E665 data fit to Eq. [| 



X 




A 


X 2 /N 


N 


0.0045(AM<7) 


0.87 ±0.16 


0.02 


2 


0.0045(^665) 




0.90 ±0.10 


0.43 


4 


0.0045(^665) 




0.94 ±0.09 


0.31 


5 


0.0080(iVA/C) 


0.75 ±0.07 


0.38 


3 


0.0080(£'665) 




0.87 ±0.10 


0.24 


4 


0.0080(^665) 




0.85 ±0.11 


1.19 


4 


0.0125(iVMC) 


0.81 ±0.05 


0.55 


5 


0.0125(5665) 


0.97 ±0.14 


1.12 


4 


0.0175(iVMC) 


0.78 ±0.06 


0.38 


5 


0.0175(£665) 


0.76 ±0.13 


0.88 


5 



"Bin center corrected from x = 0.004 
6 Bin center corrected from x = 0.005 
c Bin center corrected from x = 0.007 
d Bin center corrected from x = 0.009 
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. 3. Results from fit to CCFR data to extrapolate to 
= 0). 
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FIG. 4. F 2 (Q 2 = GeV 2 ) from different x bins. A line is 
drawn at the weighted average of all four measurements. 



TABLE II. Fit results for CCFR data. CCFR data is fit to 
Eq. 4 with A = 0.81 ±0.02 as determined by fits to NMC and 
E665 data. B, C, D, and F2 at Q 2 = results shown below. 
N = 4 for all fits. 



x B C D F V 1 {Q 2 = 0) X 2 /N 

0.0045 1.49 ±0.02 2.62 ± 0.26 0.06 ±0.17 0.04 ±0.10 OlT 

0.0080 1.63 ±0.05 2.32 ± 0.05 0.50 ± 0.05 0.22 ± 0.03 0.5 

0.0125 1.63 ± 0.05 2.39 ± 0.05 0.40 ± 0.05 0.18 ±0.03 1.0 

0.0175 1.67 ± 0.05 2.20 ± 0.05 0.65 ± 0.07 0.26 ± 0.03 0.5 
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